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Abstract 

The main obstacle in exploiting the frequency data of 5 Set stars is difficulty in mode identifica- 
tion. The 5 Set oscillation spectra, unlike those of the Sun or white dwarfs, do not exhibit very 
regular patterns. Thus, the mode identification must rely on sophisticated methods, which 
involve combined multi-passband photometry and radial velocity data, with an unavoidable 
theoretical input from stellar atmosphere models. Moreover, there are serious uncertainties in 
theory of S Set stars that have to be solved. Mode identification and determination of global 
and internal structure parameters for <5 Set stars has to be done simultaneously. I describe in 
some detail the methodology and present some recent results we obtained concerning degrees 
of excited modes, global stellar parameters, and constraints on models of subphotospheric 
convection, as well as effect of rotational mode coupling. 

Introduction 

5 Scuti stars are one of the most intensively studied group of pulsating variables. In the 
HR diagram, they are located at the intersection of the classical instability strip with the 
main sequence, and somewhat above it. It was recognised many years ago that pulsations of 
these objects, as other classical variables, are driven by the K-mechanism acting in the Hell 
ionization zone. Excited are low-order p- and g-modes with periods ranging from 0.02 d to 
0.3 d. 

Over the last 20 years, many multisite observations of these stars were carried out by 
DSN and WET networks. These campaigns have resulted in a growing number of detected 
oscillation frequencies. On the basis of these data several attempts were made towards 
construction of asteroseismic models of certain multimodal pulsators. One of such objects 
was XX Pyx, for which Pamyatnykh et al. (1998) tried to construct seismic model without the 
£ identification from photometry or spectroscopy. Another example, 9^ Tau, is a binary system 
consisting of an evolved (primary) and main sequence A-type (secondary) stars (Breger et al. 
2002), both inside the instability strip. The most multimodal and most promising object for 
asteroseismology of 5 Set stars is FG Vir. This star was studied by Guzik &d Bradley (1995), 
Viskum et al. (1998), Breger et al. (1999) and Templeton et al. (2001). Recent large 
photometric and spectroscopic campaigns, organised in the years 2002-2004 by Breger et al. 
(2005) and Zima et al. (2006), increased the number of independent oscillation frequencies 
of FG Vir to 67. In spite of all these efforts we still do not have a good seismic model for 
any 5 Set star. So far, not much has been learnt from these rich oscillation spectra. There 
are still problems with mode identification of excited modes as well as large uncertainties in 
modelling S Set pulsation to exploit the frequency data for constraining stellar model. The 
most important aspects are: turbulent convection-pulsation interaction, effects of rotation, 
mechanism of mode selection, diffusion. 

In this paper, I outline the method which gives simultaneously mode identification and 
constraints on stellar parameters and convection. I discuss also effects of uncertainties aris- 
ing from the atmospheric models and, briefly, effects of rotational mode coupling on mode 
identification. 
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Mode identification 

In the case of main sequence pulsators, the most widely used tools for mode identification are 
pulsation amplitudes and phases derived from observed variations in photometric passbands 
and in the radial velocity. If effects of rotation can be neglected, the amplitude ratio vs. phase 
difference diagrams can lead to the i degree determination, and they are independent of the 
azimuthal order, m, and the inclination angle. As was shown by Daszyriska-Daszkiewicz, 
Dziembowski &i Pamyatnykh (2003) (Paper I), in the case of 5 Set variables the photometric 
amplitudes and phases are very sensitive to the treatment of convection in the outer layers. 
This is because in calculating these observables one has to make use of the complex parameter, 
f, giving the ratio of the local flux variation to the radial displacement at the photosphere. The 
f parameter is obtained in the framework of linear nonadiabatic theory of stellar oscillation 
and, in the case of S Set stars, exhibits strong dependence on convection, as was already 
emphasised by Balona &c Evers (1999). To avoid this problem, in Paper I we invented 
a method of simultaneous determination of the £ degree and f parameter from multi-colour 
photometry and radial velocity data. The procedure consists of solving the set of observational 
equations for complex photometric amplitudes in a number of passbands, A, 
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Figure 1: Values of as a function of i for the 9.897 d ' frequency excited in /3 Cas, for M = 1.95Mq 
and three effective temperatures. In the left panel, was calculated with the Kurucz models and in the 
right panel with the NEMO2003 models. 
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Derivatives of the monochromatic flux, J^\{T^ff, \ogg), are calculated from static atmosphere 
models (Kurucz, NEMO2003, Phoenix). In general, they depend also on the metallicity 
parameter [m/H] and microturbulence velocity l;t. If the spectroscopic data exist, the above 
set of equations can be supplemented with the expression for the radial velocity (the first 
moment of line profile, M^), 

\u^R[u^ + -^v^y = Mt (2) 

In the above expressions, e is the intrinsic mode amplitude, / is the inclination angle and 
, , are disc-averaging factors weighted by the limb-darkening hx{Tgff,\ogg). For 
the limb-darkening law we use the Claret nonlinear formula. Each passband, A, yields r.h.s. 
of Eqs (1). The radial velocity data yield r.h.s. of Eq. (2). Then, the system is solved by 
the least square method assuming trial values of £. The £ identification is based on x^W 
minimisation and quantities to be determined are: £ and {ef). In Paper I we applied our 
method to three S Set stars: /3 Cas, 20 CVn and AB Cas, for their dominant frequencies. 
To this end we used amplitudes and phases in four Stromgren passbands. In all cases the 
identification of £ was unique. As an example, in Fig. 1 we plot the x^(^) dependence for 
one frequency observed in the star f3 Cas. In the two panels, the effect of using atmospheric 
models from different sources is shown. In the left panel, the x^(^) ™5 obtained adopting 
the Kurucz models, whereas in the right one, adopting the Vienna models (NE[\/lO2003). 
The method works also in the case of multiperiodic pulsators. In Daszyriska-Daszkiewicz 
et al. (2005) (Paper II) we applied the method to the most multiperiodic <5 Set star FG 
Vir. Combining the vy Stromgren photometry and radial velocity data for twelve modes, we 
arrived at a unique identification of £ in six cases, and we obtained the constraint £ < 2 m 
other six. 

The most important property of our method is that the identification of the spherical 
harmonic degree, £, is independent of any input from nonadiabatic pulsation calculations. 
Moreover, the method uses simultaneously photometry and spectroscopy by combining these 
data into the system of observational equations. For more details we refer readers to Papers 
I and II. 

Constraints on convection 

The above outlined method constitutes also a way of inferring f from observations. The 
value of f, describing bolometric flux perturbation, is determined in the pulsation driving 
zone, where the thermal time scale is comparable with the pulsation period. It means that 
this parameter is sensitive to properties of subphotospheric layers which are poorly probed 
by oscillation frequencies. In general, the f parameter depends on: mean stellar parameters, 
chemical composition, stellar convection and opacities. Thus, the strong sensitivity of the 
f parameter on convection in the case of 5 Set pulsators may be considered as an advan- 
tage. Once we know the empirical f values, we can compare them with their theoretical 
counterparts, and obtain valuable constraints on convection in subphotospheric layers. 

In Paper I, we succeeded in extracting the f parameter photometric observations for all 
studied S Set stars: f3 Cas, 20 CVn and AB Cas. We adopted Kurucz models of stellar 
atmospheres. The pulsation calculations were made assuming a simplistic approach: the 
mixing-length theory and the conveetive flux freezing approximation. In comparison of empir- 
ical f values with the theoretical ones calculated with various values of the MLT parameter 
a, we met with a problem in reproducing both the real and imaginary part of f with the same 
value of a. The general result was that observed values of were close to those calculated 
with a = 0, whereas values f/ preferred rather higher values of a. The disagreement appeared 
to be mostly correlated with the uncertainties in the atmospheric models, which I discuss in 
the next section. 
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In Paper II we applied the method of simultaneous extracting £ and f from observations 
for the most multiperiodic 5 Set star FG Vir. We relied on NEMO2003 atmosphere models. 
Combining vy Stromgren photometry and radial velocity data, we extracted the f parameter 
for twelve frequencies and compared them with the theoretical values calculated assuming two 
different treatments of convection. The first one was the standard mixing-length theory and 
the convective flux freezing approximation, as in Paper I. As the second one, we considered 
a non-local time-dependent generalisation of MLT by Gough (1977). In the first case the 
agreement was found for models with a f» 0.0, which is an evidence that convection in the 
outer layers of FG Vir is relatively inefficient. In the second case, which includes convection 
dynamics, the agreement was possible also with larger values of a, but smaller ones {a < 0.5) 
were still favoured as can be seen from Fig. 2 (taken from Paper II). 
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Figure 2: The empirical f values (dots with error bars) and the theoretical ones calculated for four values of 
the MLT parameter a, adopting a nonlocal, time-dependent formulation of MLT. The real and imaginary 
parts of f are shown in the left and the right panels, respectively. 

Modelling of S Set type pulsation with time-dependent convection treatment can be found 
also in several other papers, e.g. Grigahcene et al. 2005, Dupret et al. 2005a, Dupret et al. 
2005b. 

Uncertainties from atmospheric models 

To calculate pulsation amplitudes and phases of photometric and radial velocity variations, 
one needs an input from atmospheric models. As mentioned in the previous section, these are 
monochromatic flux derivatives over effective temperature, aj, and gravity, Qg, as well as 
the limb-darkening law, hx- In Fig. 3 we can see how non-smooth derivatives a-j-, calculated 
from Kurucz models (left panel) can produce artificial minima of derived from our method 
for a dominant mode of FG Vir. Derivatives obtained from NEMO2003 models (right panel) 
are smooth and only one minimum appears. In this case, we show also the effect of 
microturbulence velocity, on the location of minimum of x^(^eff)- The non-smooth flux 
derivatives affect also the inferred values of f. This is illustrated in Fig. 4 where empirical 
f values for 13 Cas obtained using Kurucz and Vienna models are compared with theoretical 
ones calculated for five values of the MLT parameter, a. We can see that in the case of 
Vienna models both real and imaginary parts of f are reproduced with the models assuming 
inefficient convection (q 0.0) 
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Figure 3: X a function of effective temperature for a dominant mode of FG Vir derived using Kurucz 
(left) and Vienna (right) models. Dot with error bar shows log 7"eff derived from mean colours. The right 
linked layers contain the temperature flux derivatives, Q7-, in vy Stromgren passbands. In the right panel 
the effect of microturbulence velocity is also shown. 
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Figure 4: Comparison of the f values inferred from Stromgren photometry for 13 Cas with the theoretical 
ones calculated with various MLT parameter, a. The empirical f values were obtained adopting Kurucz 
models (left panel) and Vienna models (right panel). In the right panel the effect of microturbulence 
velocity is also shown. 



Rotational mode coupling 

The most important effect of moderate rotation is mode coupling (Soufi, Goupil, Dziembowski 
1998). It takes place if the frequency difference between modes j and k is of the order of 
angular velocity of rotation, and if the spherical harmonic indexes satisfy relations: ij = ik + 2 
and mj = m^,. As elgenfunctions for Individual modes, we have to consider superpositions of 
all modes satisfying the above conditions. Hence, the photometric amplitude of a coupled 
mode Is given by (Daszyriska-Daszkiew/Icz et al. 2002) 



k 
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where the contribution from modes in non-rotating star is determined by coefficients ai, 
which are solutions from perturbation theory. Now, the location of the mode on the diagnostic 

diagrams depends on the azimuthal order, m, the inclination angle and the rotational velocity. 
We considered the stellar model with the following parameters M — 1.8Mq, log Teff = 3.866 
and log Z./Z-0 — 1.12, and the rotational velocity of about 70 km/s, which are appropriate 
for FG Vir. As an example, we consider the rotational coupling between £ — and £ — 2 
axisymmetric modes with frequencies 19.342 and 19.597 d^^, respectively. In Fig. 5 we show 
the position for coupled modes on the diagram with the Stromgren y passband and the radial 
velocity. The left panel refers to the solution dominated by the £ — component, whereas 
the right one to the solution dominated hy £ = 2. For discussion of other effects of rotation 
within the perturbative approach see e.g. Pamyatnykh (2003). 
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Figure 5: The positions of rotational coupled modes (open circles) on the Avr/ Ay vs. ^pvr — diagram. 
We considered a coupling between close £ = and 2 pair at the rotation velocity of about 70 km/s in the 
stellar model with M = 1.8Mq and log T^ff = 3.866. Filled symbols indicate positions of pure = 0, 1, 2 
modes. 



Summary 

I outlined results obtained in Paper I and II, where we proposed and applied the new method 
of simultaneous determination of the spherical harmonic degree, £, and the nonadiabatic 
parameter f from multi-colour photometry and radial velocity data. We demonstrated that 
inferring f values from such observations is possible, thus identification of £ can be done 
without a priori knowledge of f. Our method combines the photometry and spectroscopy, 
and it gives the £ identification at the highest confidence level achieved up to now. Moreover, 
by comparing empirical and theoretical f values, the method yields constraints on mean 
stellar parameters and on properties of subphotospheric layers. In the case of S Set stars, 
this is the treatment of convective transport. Inferred values of f are consistent with models 
calculated assuming rather inefFicient convection (a < 0.5). The f parameter constitutes a 
new asteroseismic tool which is complementary to oscillation frequencies. 

It is obvious that detecting more and more oscillation frequencies is of great importance, 
especially in the satellite missions time. However, it seems that asteroseismology of S Set 
stars will be served better if we focus also on those frequencies for which very accurate and 
simultaneous ground-based data from photometry and spectroscopy can be obtained. 
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